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Abstract 

One of the inevitable consequences of the ever-increasing growth of the world’s population is the 
expansion of urbanization. So, it is very important to provide a vision of the spatial development of 
cities with the aim of understanding the correct pattern of city growth and providing the necessary 
infrastructure. Since Qom metropolis faced urban growth and has recorded 95% urbanization, this 
research focused on the spatial development of urban lands around this metropolis. First, the land 
use/cover and the urban growth merit maps were produced. Land use/ land cover maps of the region 
for the years 2000, 2010 and 2020 were produced using the random forest method in the Google 
Earth Engine, and the urban growth merit map for the years 2000 and 2010 was produced using 
MCDM analyses based on GIS. Finally, based on the ANN-CA-Markov and SVM-CA-Markov 
algorithms, 2020’s land use/cover maps were simulated. The validation of the models showed that the 
SVM-CA-Markov algorithm with the ROC (0.96) was more accurate and was chosen as the optimal 
algorithm for modeling the horizon of 2040. The results indicate the increasing spatial development 
of this metropolis. The area of urban land in this region will increase from 139.62 square kilometers 
in 2020 to more than 183 square kilometers in 2040. The evaluation of the results can help the 
relevant managers in order to take necessary policies to manage the situation as best as possible. This 
importance can be realized through planning for the regular development of the road network, the 
expansion of urban green spaces, etc. official organizations and local officials should have a 
purposeful monitoring of this issue, while having full control over the development directions of this 
metropolis. 
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صص ۷۷۲ ۱۹۹ 


حكيده 


یکی از ييامدهاى اجتناب نايذير رشد روزافزون جمعيت جهان» كسترش شهرنشینی 
اسع ارايخ وى اراق متشو دارع از ترسعه فف ایی lege)‏ هدق درك الوق 
صحيح رشد شهر و فراهم آوردن زيرساختهاى لازم از اهميت بسيارى برخوردار 
است. ازآنجاکه كلانشهر قم یکی از شهرهای دركير با مسئله رشد شهرى بوده و 
آمار ٩۵‏ درصد شهرنشینی را ثبت كرده استء تمركز اين پژوهش بر واكاوى توسعه 
فضا ازاضی Gal ade‏ ابم کلان‌گهر اس برای ثبل haly giada‏ 
ورودىهاى مدل كه همان نقشههاى كاربرى/يوشش اراضى و نقشه شانستكى رشد 


شهرى منطقه هستند. توليد شدند. نقشه‌های کاربری/پوشش اراضى منطقه برای 


تاريخ دریافت: ۱۶۰۲/۰۳/۲۰ تاريخ بازنگری: ۱۸ تاريخ تصویب: ۱۶۰۲/۱۳/۰۸ 


.١‏ مقدمه 


tee‏ جغرافیا و توسعة ناحيداى شمارة چهل و چهارم 


سال‌های ۰۲۰۰۰ ۲۰۱۰ و ۲۰۲۰ با روش جنگل تصادفی در محيط سامانه Google‏ 
Earth Engine‏ و نقشه شایستگی رشد شهری منطقه برای سال‌های ۲۰۰۰ و ۲۰۱۰ به 
کمک تحلیل‌های MCDM‏ مبتنی بر GIS‏ به‌ صورت مجزا تولید شد. درنهایت اين 
نقشه‌ها وارد الگوریتم‌های ترکیبی ANN-CA-Markov‏ و SVM-CA-Markov‏ شد 
و دو نقشه برای کاربری/پوشش اراضی منطقه در سال ۲۰۲۰ شبیه‌سازی شد. 
اعتبارسنجی مدل‌ها نشان داد که الگوریتم SVM-CA-Markov‏ با مساحت زیر 
منحنی ROC‏ معادل ۰/۹۱ از صحت بیشتری برخوردار بود و برای مدل‌سازی افق 
۰ به‌عنوان الگوریتم بهینه انتخاب شد. نتایج SL‏ از توسعه فضایی روزافزون 
اين کلان‌شهر است؛ به‌طوری که وسعت اراضی شهری اين منطقه از ۱۳۹/۱۲ کیلومتر 
مربع در سال ۲۰۲۰ به بیش از ۱۸۳ کیلومترمربع در سال ۲۰۶۰ افزایش خواهد یافت. 
ارزیابی نتایج می تواند به مدیران مربوط در راستای اتخاذ سیاست‌های لازم برای 
مدیریت هرچه‌بهتر شرایط پیش رو یاری رساند. اين امر مهم می تواند از طریسق 
برنامه ریزی برای توسعه منظم شبکه معابر» گسترش فضاهای سبز شهری و.. محقق 
شود. در اين راستا سازمان‌ها و مسئولان محلی بايد ضمن اشراف کامل بر جهات 
توسعه اين کلان شهر. نظارت‌های هدفمند بر اين مسئله داشته باشند. 

کلیدواژه‌ها: تغییرات کاربری/پوشش اراضی, کلان‌شهر قم. رشد شهری. تحلیل‌های 
تصمیم گیری چندمعیاره 0160301 یادگیری ماشین» زثجیره مارکوفه سلول‌های 


خود کار. 


جمعيت جهان در گذر زمان همواره به شكل فزاینده‌ای در حال افزايش است و بهتازكى نیز از 
مرز هشت ميليارد نفر عبور كرده است (اسلام" و همکاران» ۲۰۲۳). گزارش‌ها حاكى از آن است كه 


اين روند در ادامه نيز سير صعودى خود را حفظ خواهد كرد؛ به‌طوری که پیش‌بینی می‌شود جمعیت 
حهان در سال ۰ به ۵ میلیارد نفر» در سال ۰ به W/V‏ میلیارد نفر و در سال ۰ به ۱۰/۶ 


میلیارد نفر برسد (صندوق جمعيت سازمان ملل متحد" Gis) (TY‏ اقتصادی شهرها و بهبود کیفیت 


1. Aslam 
2. UFPA (United Nations Population Fund) 


زندگی بشر سبب شده است که سهم عمده‌ای از جمعيت OL ge‏ در شهرها سکنا گزینند (گوپتا و 
بهارات ‏ ۲۰۲۲). در حال حاضر» 00 درصد جمعیت Ob ge‏ در شهرها زندگی می‌کنند و پیش‌بینی 
می‌شود که همسو با افزايش جمعیت جهان. اين رقم نیز افزايش يابد و در سال ۲۰۵۰ به ۷۰ درصد 
برسد (صندوق جمعیت سازمان ملل متحد. ۲۰۲۲). افزايش روزافزون جمعیت شهرنشین و به دنبال 
آن رشد شهرها؛ منجر به پیشرفت جوامع بشری و تقویت اقتصاد ملی کشورها خواهد شد 
(موسوینگان" و همکاران. ۲۰۲۰). اما توسعه فضایی شهر می‌تواند پیامدهای بی‌شماری ater jl‏ 
افزايش دمای سطح زمین (LST‏ (گوپتا؛ و همکاران. ۲۰۲۳). تخریب و از دست دادن تنوع زیستی 
(مک‌دونالد" و همکاران. ۲۰۲۰) تراکم جمعیت و حاشیه‌نشینی (پیله ورأء ۲۰۲۰ و تغییرات در 
کاربری/پوشش اراضی LULCY‏ را به همراه آورد (گوان" و همکاران» ۲۰۲۲). 

تغييرات LULC‏ نحوه تعاملات بين انسان و محيط پیرامون Ol‏ را در طول زمان نشان می‌دهد 
(مظفری‌پور و همکاران. ۲۰۲۲). اين پدیده اثرات گسترده‌ای در مقیاس‌های محلی, منطقه‌ای و جهانی 
دارد (رگسا و نونز“ ۲۰۲۲). تهدید امنيت غذایی (لونیولو"" و همکاران ۲۰۲۱ کاهش كيفيت 
خدمات اکوسیستم (فنگ '' و همکاران. ۲۰۲۲) و بیابان‌زایی (مینائی و همکاران. ۲۰۱۸ نمونه‌ای از 
پیامدهای زیان‌بار تغییرات LULC‏ هستند. با توجه به اينكه از یک سو تغيير الگوی کاربری زمین 
اثرات گسترده‌ای بر زیستگاه‌های طبیعی دارد (حسین و همکاران. ۲۰۲۳) و از سوی دیگر بررسی 
سناریوهای متحمل بر آینده تغييرات LULC‏ حاکی از گسترش فزاینده اراضی شهری است 
“Esl +)‏ و همکاران» ۲۰۱۹ ارزیابی روند تغییرات LULC‏ یکی از اصلی‌ترین اولويتهاى 
برنامه‌ریزان محیط‌زیست و متولیان امور شهری محسوب می‌شود (گیبس و مایورک"؛ ۲۰۲۰)؛ 


ازاین‌رو در دهه‌های اخير توسعه فضایی ناموزون اراضی شهری. به یکی از موضوعات جهانی تبدیل 


. Gupta & Bharat 

. Musavengane 

. Land Surface Temperature 
. Gupta 
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. Land Use Land Cover 
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شده و توجه مدیران و برنامه‌ریزان شهری در کشورهای توسعه‌یافته و درحال‌توسعه را به خود جلب 
کرده است (هارماى' و همکاران» ۲۰۲۱). در همین راستاء يايش جهات توسعه شهر و پیش‌بینی الگوی 
توسعه فضایی شهرها در آینده که برای تحقق اهداف توسعه پایدار نيز ضروری است. مدنظر متولیان 
حوزه شهری قرار گرفته است (ما" و همکاران ۲۰۲۲). 

گفتنی استء اين موضوع در کلان‌شهرها که سهم عمده‌ای از رشد شهری یک کشور را به خود 
اختصاص مىدهندء از اهمیت دوچندانی برخوردار است (ابراهیم و همکاران» ۲۰۲۱). کلان‌شهر قم 
یکی از کلان‌شهرهای رسمی ايران است که در ساليان اخير جمعيت افراد شهرنشين در OF‏ به‌شدت رو 
به افزايش بوده است. در آینده نيز می‌توان ادامه همین روند را برای اين کلان‌شهر متصور شد (مرکز 
آمار col yl‏ ۱۳۹۵). بدیهی است که اين امر سبب توسعه فضایی بیش‌ازپیش اين کلان‌شهر خواهد شد. 
بدون تردید معرفی روشی برای مدل‌سازی الگوهای مکانی-زمانی تغييرات کاربری و پوشش آینده 
اراضی در اين کلان‌شهر بسیار ضروری به نظر می‌رسد؛ ازاین‌رو مطالعه حاضر سعی دارد با ارائه 
رویکردی نوین» به مدل‌سازی توسعه فضایی کلان‌شهر قم در افق ۲۰۶۰ بپردازد و مدیران مربوط را 
قادر سازد تا علاوه‌بر جلوگیری از جه كيرى فضایی ناموزون اين کلان‌شهر زیرساخت‌ها و 


تمهیدات لازم برای گسترش فیزیکی آن را نيز فراهم آورند. 


مدل‌سازی وسعه فضایی شهرها به‌عنوان یکی از اصلی‌ترین چالش‌های مطالمات شهری از spe‏ 
توجه پژوهشگران را به خود جلب کرده است (كيم' و همکاران ۲۰۲۲). در اين ميان امروزه 
بهره‌مندی از علم سنجش از دور“ و سیستم اطلاعات جغرافیایی (GISY‏ و همچنین الگوریتم‌های 
یادگیری ماشین (ML‏ در زمره روش‌های کارآمد برای مدل‌سازی رشد فضایی شهرها به شمار 
می‌رود (هایند" و همکاران. ۲۰۲۲). یکی از مدل‌های محبوب در زمینه مدل‌سازی توسعه فضایی 
شهرها که انطباق‌پذیری بسیار زیادی با سنجش از دور و GIS‏ دارده مدل ترکیبی CA®Markov‏ 


. Harmay 

. Ma 

Kim 

. Remote Sensing 

. Geographic Information System 
. Machine Learning 

Hind 

. Cellular Automata 


س وح ين خي بي ي نف ص 


اين مدل روند توسعه فضایی شهرها را شبیه‌سازی کرده‌اند؛ به عنوان نمونه‌ای از اين مطالعات می‌توان 
به پژوهش (سنگ" و همکاران ۰1۱( در ناحيه شهری فانگشان " در يايتخت کشور جين اشاره کرد 
كه در آن ابتدا به کمک داده‌های ماهواره‌ای سنجش از دور نقشه چند زمانه LULC‏ منطقه تولید 
شده است. سپس در ادامه آن‌ها با استفاده از مدل ترکیبی CA-Markov‏ نقشه LULC‏ منطقه را 
برای سال ۲۰۱۵ شبیه‌سازی کردند. نتایج حاصل از شبیه‌سازی آن‌ها حاکی از گسترش بیش‌ازپیش 
مبتنی بر سنجش از دور و GIS‏ روند تغییرات LULC‏ مناطق شهری دره كاتماندو” واقع در کشور 
نيال را بين سال‌های ۱۹۸۸ تا ۲۰۱۱ بررسی کردند. شواهد نشان داد که در اثنای اين مدت» وسعت 
مناطق شهری اين محدوده با رشد سالانه ۹/۱۵ درصد از ۳ كيلومترمربع به ۱/۳۵ کیلومترمربع 
افزايش يافته است. درعين حال شبیه‌سازی تغییرات LULC‏ منطقه برای سال‌های ۲۰۲۶ و ۲۰۳۲ که با 
استفاده از _ CA-Markov‏ انجام شده است نيز ادامه همین روند را پیش‌بینی می‌کند؛ به‌طوری که 
نتایج اين پژوهش OL‏ داد كه وسعت اراضی شهری در سال‌های مذکور به‌ترتیب به ۲۰۰ و ۲۳۸ 
کیلومترمربم خواهد رسید. نقشه LULC‏ شبیه‌سازی‌شده برای سال ۲ همجنين گویای اين امر بود 
که در ازای اين میزان افزایش اراضی شهری. حدود ۵۰۰ کیلومترمربع از اراضی کشاورزی تخریب 
خواهد شد. پژوهش کیسابا و لى' (۲۰۲۳) در شهر دودما" در کشور تانزانیا ازجمله مطالعات اخير در 
زمينه مدل‌سازی توسعه فضايى شهرها به حساب می آید که با روش CA-Markov‏ انجام شده است. 
نتایج حاصل از شبیه‌سازی تغييرات LULC‏ اين شهر برای سال ۲۰۳۰ نشان می‌دهد كه مساحت 
اراضی شهری در اين منطقه به میزان ۳۹/۹۲ درصد افزایش خواهد یافت. همسو با اين مطالعات. 
شفیع زاده‌مقدم و “eels‏ (۰۲۰۱۳ رهنما (YYA)‏ قا“ و همکاران (۲۰۲۱) و طارق ٠١‏ و همکاران 


AVY)‏ نيز در پژوهش‌های خود به کمک مدل CS sa CA-Markov‏ به مدل‌سازی توسعه 


1. Fu 
2. Sang 
3. Fangshan 
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6. Kisamba & Li 
7. Dodoma 

8. Shafizadeh Moghadam & Helbich 
9. Baqa 

10. Tariq 


alee €‏ جغرافیا و توسعة ناحیهای شمار؛ چهل و چهارم 


فضایی شهرهای بمبئى ' کراچی " و پیشاور" يرداختند. مدل‌سازی‌های حاصل از اين مطالعات همكى از 
گسترش دوچندان اراضی شهری خبر می‌دهند. 

در ايران نيز مطالعات متعددی در زمینه مدل‌سازی توسعه فضایی cle gt‏ مختلف با تکیه بر مدل 
CA-Markov‏ انجام شده است كه از جمله آن‌ها می‌توان به پژوهش امینی پارسا و همکاران 
GYAL)‏ اشاره کرد که در آن با استفاده از این مدل و نقشه‌های LULC‏ منطقه که برای سه دوره 
۷ ۰ و :۲۰۱۶ agi‏ شده بود. کاربری اراضی شهرستان ملکان در استان آذربایجان غربی برای 
٤‏ سال آينده مدل‌سازی شده است. نتایج اين تحقیق نشان داد که در دوره زمانی ۲۰۱۶ تا ۲۰۲۸ 
کاربری‌های مسكونى و کشاورزی» كاربرى غالب در اين منطقه خواهند بود. در نمونهاى دیگر» میثاق 
و همكاران ATAD‏ با تلفيق تحليلهاى تصميم كيرى جندمعياره مبتنى بر GIS‏ با مدل ترکیبی CA-‏ 
Markov‏ روند گسترش افقى کلان‌شهر تبريز را برای سال ١407‏ مدل‌سازی كردند. نتایج حاصل از 
مدل‌سازی Lol‏ نشان‌دهنده رشد ۲۵ درصدی اراضی شهری اين کلان‌شهر بود. همسو با اين مطالعات؛ 
اکبری و همکاران (۱۳۹۸) نيز پژوهشی مشابه داشتند و توسعه فضایی شهر رشت را برای سال ۲۰۲۷ 
شبیه‌سازی کردند. نتایج حاصل از تحقیق آن‌ها مؤيد افزايش مضاعف وسعت اراضی شهری در این 
منطقه بود. در همین سال» یوسف‌زاده و همکاران نیز از روش مذکور برای بصری‌سازی توسعه فضایی 
کلان‌شهر مشهد و حومه تا سال ۲۰۵۰ استفاده کردند. نقشه شبیه‌سازی‌شده نشان داد که گسترش 
فضايى كلانشهر مشهد عمدتاً به سمت غرب خواهد بود. از طرفی در اين منطقه رشد مناطق 
انسانساخت با كاهش پوشش گیاهی همراه خواهد شد. 
ازآنجاکه پدیده رشد شهری فرایندی غیرخطی است و الگوریتم‌های ML‏ در تشخیص روابط غير 
خطی بسیار کارآمد هستند (کیم و همکاران, ۲۰۲۲ اخيراً اين الگوریتم‌ها در كنار مدل CA-‏ 
۷ و در ترکیب با مدل CA‏ در فرایند مدل‌سازی توسعه فضایی شهرها استفاده شده‌اند. از 
ميان تحقیقاتی که با استفاده از الگوریتم‌های یادگیری ماشین انجام شده است. می‌توان به پژوهش 
القاضی * و همکاران (۲۰۲۱) اشاره کرد که در آن پس از طبقه‌بندی نقشه‌های چند زمانه LULC‏ به 


کمک الگویی تلفیقی از الگوریتم یادگیری ماشین شبکه عصبی مصنوعی و مدل Glad gh‏ خودکار 


1. Mumbai 
2. Karachi 
3. Peshawar 
4. AlQadhi 


سال بيست و يكم بهبود مدل‌سازی توسعه فضایی شهرها با تلفیق ... ۳9 


(ANNI-CA)‏ توسعه فضايى شهر ابها" در عربستان را برای سال ۲۰۲۸ شبیه‌سازی کردند. نتايج 
شبیه‌سازی شواهدى مبنى بر رشد بىرويه وسعت اراضى شهرى در اين منطقه ارائه كرد. در پژوهشی 
مشابه با مطالعه قبلی. خان و سودهیر " (۲۰۲۲) برای يايش توسعه فضايى شهر اسلام‌آباد در كشور 
پاکستان از همین مدل (ANN-CA)‏ بهره جستند. شباهت فضایی VO‏ درصدی نقشه حاصل از این 
مدل با واقعیت زمینی» سبب شد که آن‌ها از اين رویکرد برای شبیه‌سازی روند تغییرات کاربری 
اراضی منطقه برای سال‌های ۲۰۳۱ و ۲۰۶۱ نيز استفاده کنند. نتایج اين مطالعه حاکی از تخریب 
جنگل‌های منطقه و تبدیل Of‏ به اراضی شهری است. پژوهش شفیع‌زاده مقدم و همکاران (۲۰۲۱) 
نمونه‌ای از مطالعات اين چنینی است که در ايران انجام شده است. در اين مطالعه» نویسندگان به 
کمک نقشه‌های کاربری/پوشش زمين در سال‌های ۰۱۹۹۸ ۲۰۰۸ و ۲۰۱۸ و ترکیب الگوریتم انتخاب 
ویژگی رو به جلو پیش‌رونده با الگوریتم جنگل تصادفی * (FFS-RF)‏ و همچنین تلفیق آن با 
مدل CA‏ نقشه تغییرات LULC‏ منطقه را برای بازه زمانی بين سال‌های ۱۹۹۸-۲۰۰۸ و ۲۰۱۸- 
۸ تهيه كردند. شبيهسازى تغييرات LULC‏ اين محدوده OLE‏ داد که تا سال ۲۰۰ تهران. کرج و 
شهرهاى مجاور آن از نظر فضايى به هم متصل می‌شوند و ابرشهر غول ييكر* را تشكيل خواهند داد. 
بررسى بيشينه تحقيق نشان مىدهد كه بيشتر مطالعات در زمينه موردبحثء با روش‌های مبتنى بر 
سنجش از دور و GIS‏ و همجنين بهره‌مندی از مدل تركيبى CA-Markov‏ يا تركيب یکی از 
الكوريتمهاى ML‏ با مدل CA‏ انجام شده است. اما اين مطالعه سعى دارد كه با نكاهى نوآورانه و با 
هدف بهبود مدل‌سازی توسعه فضايى شهرهاء اين فرايند را با تلفيقى از اين دو روشء يعنى با 
الگوریتمی تركيبى مبتنى بر ML-CA-Markov‏ انجام دهد. بهره‌مندی از روش BWM‏ در زمينه 
موردبحث و تهيه نقشه چندزمانه LULC‏ منطقه در سامانه متن باز Google Earth (GEE)‏ 
Engine‏ را كه براى پردازش تصاوير ماهوارهاى به سختافزارهاى يبشرفته و كرانقيمت نياز ندارد و 
تصاوير تصحیح شده رابه صورت بازتاب زمينى (Surface Reflectance)‏ در اختيار كاربران قرار 
می‌دهد. می‌توان از So‏ ابتکارات نوآورانه اين تحقیق برشمرد. شایان ذکر است که پژوهش حاضر 


اولين تحقیقی است كه به صورت جامع تغییرات آتی LULC‏ کلان‌شهر قم را واکاوی کرده است. 


1. Artificial Neural Network 
2. Abha 

3. Khan & Sudheer 

4. Random Forest 

5. Gigantic City 


alee ea‏ جغرافیا و توسعة ناحیهای شمارة چهل و چهارم 


۳ روش شناسى تحقیق 

در اين پژوهش به‌منظور مدل‌سازی توسعه فضایی کلان‌شهر قم ابتدا از طریق مطالعات اسنادی و 
کتابخانه‌ای, عوامل محرک توسعه فضایی شهر در منطقه موردمطالعه استخراج شدند. در ادامه 
نقشه‌های دیجیتالی هرکدام از اين معیارها در محیط GIS‏ تهیه شده و سپس بابه‌کارگیری روش 
71 که یکی از روش‌های نوين MCDM‏ مبتنی GIS‏ است (مینائی و همکاران. ۲۰۲۱ هرکدام 
از معیارهای انتخابی وزن‌دهی شده و درنهایت نقشه شایستگی رشد شهری منطقه مورد مطالعه تولید 
شده است. یکی دیگر از ورودی‌های لازم برای انجام اين پژوهش. نقشه‌های چندزمانه LULC‏ 
هستند که در پژوهش حاضر اين نقشه‌ها با الگوریتم RF‏ و در سامانه Google Earth Engine‏ 
تولید شده‌اند. در ادامه يس از تهيه ورودی‌های مدلء یک‌بار با الگوریتم ترکیبی ANN-CA-‏ 
۷ و بار دیگر با ترکیب SVM-CA-Markov‏ تغییرات LULC‏ منطقه برای سال ۲۰۲۰ 
شبیه‌سازی خواهد شد. گفتنی است که الگوریتم‌های یادشده در محيط نرم افزار 1617566 پیاده‌سازی 
شده است. در انتها صحت‌سنجی شبیه‌سازی با استفاده از منحنی مشخصه عملکرد نسبی ) CROC!‏ 
انجام شده و الگوریتم بهینه برای شبیه‌سازی LULC‏ کلان‌شهر قم در سال ۲۰۶۰ انتخاب خواهد شد. 


فرایند انجام اين پژوهش در شکل ۱ آمده است. 


1. Receiver Operating Characteristic 


بال سك ويك بهبود مدلسازى توسعه فضايى شهرها با تلفيق ... Y.V‏ 


مطالعه اسنادی پیرامون فاکتورهای مؤثر | 


بر رشد شهری 


نقشه شبیه‌سازی شده رشد نقشه شبیه‌سازی شده رشد 
شهری سال۲۰۳۰ با شهری سال۲۰۳۰ با 
ANN-CA-Markoy SVM-CA-Markovy‏ 


نقشه شبیه‌سازی شده رشد شهری 
برای سال ۲۰۴۰ 


شکل ۱. فرایند روش شناسى پژوهش 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


.١ 1‏ منطقه موردمطالعه 


محدوده مطالعاتى اين پژوهش دربرگیرنده حريمى در پیرامون كلانشهر قم است که با مساحتى 
بالغ بر ۱۰3۵ كيلومترمربع در مركز كشور ايران واقع شده است (شكل ۲). تويوكرافى اين محدوده 
شامل يهنهاى بسيار كم شيب با حداكثر ارتفاع 11٠١‏ متر است كه در مخروطه افكنه قمرود قرار دارد 
(رنجبر سعادت آبادى و فتاحی» ۱۳۹۷ ص. 1۳). آب و هواى منطقه نيز نيمهبيابانى با تابستان‌های گرم 
و خشک و زمستان‌های کم و بيش سرد است. از طرفى ميانكين سالانه بارش و دماى اين منطقه به 
ترتيب ۱۳۸/۹میلی‌متر و ۱۸/6 درجه سلسيوس ثبت شده است (اداره كل هواشناسى استان قم VEY‏ 
از نظر ذخاير معدنى نيز اين منطقه. درب ركيرنده بزرگ‌ترین معدن منگنز خاورميانه (معدن ونارج) و 
معادن متعدد سرب» روی, آهن و مس است (رفعتى و همکاران ۱۳۹۶ ص ۱۵۵). گفتنی است که 
قسمت عمده اراضی شهری اين منطقه را کلان‌شهر قم تشکیل می‌دهد. اين کلان‌شهر که هشتمین 
کلان‌شهر رسمی ايران شناخته می‌شود. طی چند دهه اخیر تغییرات جمعیتی چشمگیری را تجربه کرده 
است و جمعیت آن از ۹۹4۹٩‏ نفر در سال ۱۳۳۵ به ۱۲۰۱۱۵۸ نفر در سال ۱۳۹۵ افزايش يافته است 


(مرکز آمار ايران» ۱۳۹۵). گسترش مهاجرت را می‌توان عامل اصلی رشد جمعیت در اين منطقه 


ales Y-A‏ جغرافیا و توسعة ناحیه‌ای شمارة چهل و چهارم 


دانست (مجتبایی و نوابخش» ATMA‏ ص ۲۱۱). گردشگری مذهبی و تسهیل راه‌های مواصلاتی به اين 
کلان‌شهر را که نقش گذرگاهی در حاشیه کویر مرکزی ايران دارد. می‌توان از دیگر عوامل مؤثر در 
رشد بی‌رویه جمعیت در اين منطقه برشمرد (رفعیان و زاهد. ۹ ص. AVY‏ پیامد اصلی جنين 
وضعیتی» رشد سریع کالبد شهر و گسترش شهرنشینی است؛ به‌طوری که هم‌اکنون نسبت شهرنشینی در 
اين کلان‌شهر بیش از ۹۵ درصد است (مرکز آمار ایران» ۱۳۹۵). 
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شكل ۲. منطقه موردمطالعه 


مأخذ: ترسیم براساس یافته‌های پژوهش ۱۶۰۲ 


۳ ۲ عوامل محرک رشد شهری 
در این پژوهش شنت معیار فاصله از حاده‌هاء فاصله از اراضی بای فاصله از اراضی زراعی» 
فاصله از مناطق انسان‌سانحت, فاصله از گسل» شیب ارتفاع و فاصله از مناطق صنعتی از طریق مطالعه 


ا کت وگ بهبود مدل‌سازی توسعه فضایی شهرها با Guild‏ ... ۲۹ 


جدول ۱. عوامل محرک رشد شهری 
مأخذ: اسناد مربوطه 


فأضلة از Wale‏ (كيو' و همکاران ۲۰۲۲) فاصله از كسل 
فاصله از اراضى باير (الشریفی ۲ و همکاران ۲۰۲۲) (مارکوس E‏ و همکاران ۲۰۲۱) 


۱ (حسین شبهو و اسلام ؛ ۱ 
فاصله از مناطق انسان‌ساخت فاصله از مناطق صنعتی (زنگ" و همکاران» ۲۰۲۲) 
(oY‏ 


در ادامه با توجه به ماهيت متفاوت هرکدام از اين معیارها و نظر به اینکه روش BWM‏ صرفاً 


ی وای اة فاي اه فا شر ها etal‏ ااا E‏ كدان lo‏ ت 
هر اد به ا ادو کو ال انردق و سهان بک و OR‏ ني 
صفر و ۱ تعریف می‌شود. در اينجا حداكثر مقادير برای معيار فاصله از گسل مطلوب است که 
به‌صورت زیر بیان می‌شود: 

× i = (xij) / Oj - xj) O) رابطه‎ 


x! ij کت‎ (xij => at) / og = xj™) (Y) رابطه‎ 


1. Qiu 

2. Kamranzad 

3. Alsharif 

4. Marques 

5. Jana 

6. Ding 

7. Hossain Shubho & Islam 
8. Zhang 


alee 11۳‏ جغرافیا و توسعة ناحیهای شمارة چهل و چهارم 


شکل ۳ لایه‌های استانداردشده استفاده‌شده در فرایند مدل‌سازی رشد شهری سال ۲۰۲۰ 
( الف: فاصله جادههاء ب: فاصله از اراضی بای پ: فاصله از اراضی زراعی» ت: فاصله از اراضی انسان‌ساخت. 
ث: فاصله از گسل» ج: شیب چ: ارتفاع» ح: فاصله از مناطق صنعتی) 
مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


۳ ۳. نقشه‌های کاربری/پوشش اراضی 

برای ag‏ نقشه‌های LULC‏ کلان‌شهر قم از تصاویر ماهواره Landsat‏ در سال‌های ۲۰۰۰ 
۰ و ۲۰۲۰ در ماه‌های بدون ابر و با بیشترین میزان سبزینگی استفاده شده است (جدول AY‏ 
گفتنی است که در پژوهش حاضرء تصاوير دوره زمانی ۲۰۰۰ تا ۲۰۱۰ برای کالیبراسیون مدل و 
تصاویر دوره زمانی ۲۰۱۰ تا ۲۰۲۰ برای فرابند شبیه‌سازی در سال ۲۰۶۰ استفاده شده است. در اين 
راستاء ابتدا تصاویر مذکور در محيط سامانه GEE‏ فراخوانی شده و سپس به کمک الگوریتم 
طبقه‌بندی نظارت‌شده جنگل تصادفی (بریمن ؛ ۲۰۰۱ با وارد كردن ۳۰۰ نقطه تعلیمی" نقشه‌های 
LULC‏ منطقه در ۵ كلاس زراعت و باغات» Sol‏ جنگلی» پهنه‌های آبی. انسان‌ساخت و ple‏ برای 


سالهاى ياد شده تولید peer rye oe‏ (شکل nes‏ 


1. Breiman 
2. Training Point 


ارزیابی صحت اين نقشه‌ها با تعداد ۰ نمونه تست" و دو سنجه‌ی صحت کلی" و ضریب كايا" 
انجام شد. صحت JS‏ طبقه‌بندی برای تصاویر مربوط به سال‌های و ۰ به‌ترتیب 
1 درصدء ۹۵/۵7 درصد و AVIE‏ درصد به دست آمده است. ضريب كايا نيز برای سال‌های 
مذکور به ترتیب برابر با ۰/۸۹ ۰/۸۸ و ۰/۹۰ است. اطلاعات مربوط به هرکدام ازاين نقشه‌ها در 


جدول ۲ آمده است. 


جدول ۲. ویژگی‌های تصاوير ماهواره‌ای Landsat‏ 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


شماره ماهواره سنحنده 1 قدرت تفکیک مکانی (متر) 
تصوير 1 
Yr ANN‏ 
ETM Landsat \‏ لى 
7 ۲۰/۷/۰۱ 
1/7/۰۱ 
ETM Landsat y‏ لى 
7 ۲1/۷/۰۱ 
L nd‏ ۲۰/7/۰۱ 
OLI andsat 3‏ 
8 ۲۲/۷/۰۱ 


کاربری | پوشش اراضی سال 


or‏ وكا 
1 وكا 
1 وكا 


س زراعت و باغات 9 

فق يارك جنگلی 3 

1 يهندهاى آبى‎ ml 

pa انسان ساخت‎ m 

3 2: *9 
0 4 8 16 24 

كيلومتر 


شكل £ نقشه‌های كاربرى / يوشش اراضى توليد شده در سامانه GEE‏ 


مأخذ: یافته‌های يؤوهشء ۱۶۰۲ 


1. Testing Point 
2. Overall Accuracy 
3. Kappa 
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مأحذ: یافته‌های پژوهش. ۱۶۰۲ 
کاربری / پوشش 


"۳ کاربری / یوشش اراضی ۱۳۷۹ کاربری / پوشش اراضی ۱۳۸۹ کاربری / پوشش اراضی ۱۳۹۹ 
راصى 


LULC 


| ۲۵۳ | Y 17/14 | ۲۵۸ | nm | we | wAY | زراعت | زراعت و باغات‎ 


د هه و Se‏ د 


۳ ۶ روش BWM‏ 
یکی از روش‌های MCDM‏ مبتنی بر GIS‏ که به‌تازگی مدنظر محققان قرار گرفته است. روش 
بسیار ساده‌تر است و مقایسات زوجی کمتری دارد که اين امر سبب می‌شود وزن‌دهی قابل اعتمادتر 
شود (مینائی و همکاران» ۲۰۲۱). در اين روش از دو بردار مقایسه زوجی برای تعیین وزن معیارها 
استفاده می‌شود. ابتدا بهترين (مطلوب ترين» مهمترين) و بدترین (نامطلوب‌ترین. کم‌اهمیت‌ترین) معیار 
توسط تصمیم كير نده شناسایی شد و يس از آن» مقایسه‌های زوجی بين هر یک از این دو معیار 
شناسایی وزن‌ها استفاده شده است (مينائى و همکاران» ا مراحل روش BWM‏ که در Lew!‏ 

برای وزد‌دهی عوامل محرک رشد شهرى استفاده شده است» به شرح زیر است: 
مرحله ۱. شناسایی و استخراج معیارهای تصمیم‌گیری C2, ..., Cn}‏ ,61)؛ 
مرحله ۲. انتخاب بهترین و بدترین معیار از ميان معیارهای موجود. در اين مرحله. تصمیم‌گیرنده 


بهترین و بدترین معیار را به‌صورت کلی مشخص می کند؛ 


1. Best-Worst Multi-Criteria Decision-Making Method (BWM) 


مرحله ۳. مقایسه زوجی معیارها. بدین‌منظور» تصمیم‌گیرنده میزان ارجحیت بهترین معیار را 
نسبت به plo‏ معیارها با اعداد ۱ تا ٩‏ مشخص می کند که در OT‏ عدد ٩‏ نشان‌دهنده بیشترین ارجحیت 
و عدد ۱ نشان‌دهنده کمترین ارجحیت است. نتيجه اين مقایسه برداری به نام Best-to-others‏ 
۲ را تشکیل می‌دهد که بەصورت AB = (481, 282, ..., ABn)‏ نمایش داده می‌شود. در اين 
بردار» 28 نشان‌دهنده ارجحیت بهترین معیار (B)‏ نسبت به معیار (J)‏ است که BBEL‏ 

مرحله ۶. در اين مرحله» تصمیم‌گیرنده همان فرایند مرحله سوم را برای بدترین معیار انجام 
می‌دهد. نتيجه مقایسه زوجی معیارها در اين مرحله برداری به نام Other-to-worst vector‏ را 
تشکیل می‌دهد که به صورت Aw = (aiw, 82۷ ..., anw)!‏ نمايش داده می‌شود. در اين بردار؛ 
ajw‏ نشان‌دهنده ارجحیت معیار ([) نسبت به بدترین معیار (W)‏ است که ۵«<1؛ 

مرحله ۵. وزن‌دهی بهینه معیارها Wn’)‏ و (Wr, W2",‏ در این مرحله به‌منظور تعیین وزن 
بهینه هرکدام از daglas‏ زوج‌های a = apy‏ و diw‏ = = تشکیل می‌شود. سپس برای برآورده كردن 
اين شرایط در همه معيارهاء بايد راه‌حلی بيدا شود تا عبارت an|‏ = و [E-a‏ را برای همه 
معیارهایی که حداقل شده است. حداکثر کند. با توجه به مثبت بودن اوزان معیارها می‌توان مدل را 
به‌صورت رابطه (۳) فرموله کرد: 


: wa "i 
min max; ۳ -al.l = aw )۳( رابطه‎ 


1 


Sil 


wj = «Û for all j 
همچنین می‌توان مدل بالا را به‌صورت رابطه (4) نيز بیان کرد.‎ 
miné (£) رابطه‎ 
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پس از به دست آوردن مقادیر بهینه» بايد نرخ سازگاری در روش BWM‏ محاسبه شود. نسبت 
سازگاری BWM‏ با استفاده از GF‏ و رابطه )0( به دست می‌آید (مینائی و همکاران» ۲۰۲۱). مقدار *ي 
بيشتر نشاندهنده نرخ سازگاری بيشتر است. شاخص سازكارى در جدول (E)‏ ام است. 

رابطه )0( = 


Consistency Ratio = Caia ia 


جدول؛. شاخص سازگاری در روش BWM‏ 


alee 


شاخص 
mega 0 Soc meee‏ 

سازگا 

زداری 


۳ ۵. شبکه عصبی مصنوعی 
یکی از الگوریتم‌های ML‏ مرسوم در مدل‌سازی توسعه فضایی شهرها؛ الگوریتم شبکه عصبی 
مصنوعی (ANN)‏ است (چتری و سوراورا. ۲۰۲۱). روش شبکه عصبی پرسپترون چندلایه 
(MLP)‏ که روملهارت" و همكاران OAN‏ ارائه كردند. یکی از يركاربردترين شبكههاى عصبى 
مصنوعى شناخته می‌شود. اين الكوريتم شامل يك لايه ورودی» یک لايه خروجى و یک يا چند لايه 
پنهان* است (دانشفر و همكاران. ۲۰۲۰). نورون‌ها عنصر اصلی تشکیل‌دهنده لايهها در اين شبكه 
پیش‌رونده هستند (ياريان* و همکاران» ۲۰۲۲). اوزان در شبكه عصبى مصنوعى با به‌کارگیری توابع 
فعال‌سازی به نورون‌های خروجى منتقل می‌شوند (گریما" و همکاران» ۲۰۲۲). خروجىهاى حاصل از 
هرکدام از اين توابع در بازه‌های خاصى قرار می‌گیرند. در پژوهش پیش‌رو در معماری الگوریم 
ANN‏ از یک تابع فعال‌سازی ٩12۳0010‏ (هرینگتون ۱۹۹۳) که برد آن بين بازه صفر تا یک است و 
تعداد چهار a‏ پنهان به‌منظور تولید نقشه پتانسیل تبدیل کاربری/پوشش اراضی LULC)‏ استفاده 


awl odd 


1. Chettry & Surawar 

2. Multilayer Perceptron Neural Network 
3. Rumelhart 

4. Hidden Layer 

5. Yariyan 

6. Girma 


VY‏ ماشین بردار پشتیبان 


در اين پژوهش برای دستیابی به الگوریتم بهینه برای مدل‌سازی توسعه فضایی کلان‌شهر قم 
علاوه‌بر الگوریتم ANN‏ از الگوریتم SVM‏ که کورتس و واپنیک" )1990( معرفی کردند نیز در 
قسمت تولید نقشه پتانسیل تبدیل کاربری/پوشش اراضی LULC)‏ استفاده شده است. اين الگوریتم 
مبتنی بر هوش مصنوعی ازجمله روش‌های طبقه‌بندی نظارت‌شده است که به کمک مرز تصمیم گیری 
(ابر صفحه) " کلاس‌ها را از یکدیگر تفکیک می‌کند (كيم و همکاران. ۲۰۲۲). در اين الگوریتم. توازن 
ميان حاشیه " و حطای طبقه‌بندی توسط عدد ثابت © کنترل می‌شود (باجراچاریا و سلطاناگ ۲۰۲۲).اگر 
طبقه‌بندی داده‌ها با استفاده از یک تابع خطى امکان‌پذیر نباشد. ماشین بردار پشتیبان می‌تواند با استفاده 
از توابع کرنل» به‌عنوان طبقه‌بندی‌کننده غیرخطی استفاده شود. از طرفی انواع مختلفی از توابع کرنل 
وجود دارد اما تابع پایه شعاعی (RBE)‏ به‌عنوان بهترین انتخاب از بين دیگر توابع کرنل گزارش 
شده است (فنگ" و همکاران ۲۰۱۳). به‌طور کلی در SVM‏ يس از طبقه‌بندی می‌توان به کمک یک 
تابع RBF‏ پتانسیل تبدیل تغییرات LULC‏ و توسعه مناطق شهری را تخمین زد (هلدر" و همکاران, 
۳ كه در اینجا نیز از همین تابع برای تولید نقشه‌پتانسیل تبدیل تغییرات LULC‏ استفاده شده 


است: 
۳ ۷ مدل CA-Markov‏ 


در اين پژوهش برای مدل‌سازی تغییرات آتی LULC‏ کلان‌شهر قم و حریم پیرامون آن از مدل 
تلفیقی CA-Markov‏ استفاده شده است. مدل Markov‏ زنجیره‌ای دنباله‌دار از فرایندهای تصادفی 


است که حالت‌های تبدیل و نرخ تبدیل بين انواع کاربری‌ها را به‌صورت کمی ارائه می‌دهد. خروجی 
مدل Markov‏ ماتریسی mxm‏ از احتمال تغییرات در بازه زمانی مدنظر است که می‌توان Ol‏ را 


1. Cortes & Vapnik 

2. Hyperplane 

3. Margin 

4. Bajracharya & Sultana 
5. Radial Basis Function 
6. Feng 

7. Halder 
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Bi Fx Po 
Leey - ير‎ ele ی و‎ < Po Pe یگ‎ CV) رابطه‎ 
Paa Pow Fan 


در رابطه بالا رو و بمب به ترتیب حالت‌های کاربری زمین در زمان tH)‏ و ا است. 


ازانجاکه در خروجی مدل Markov‏ هیچ مولفه مکانی وجود ندارد ( ایسین كارالار' و همکاران؛ 
۲ ) ضرورى است که برای تکمیل فرایند مدل‌سازی از یک مدل مکمل مانند مدل CA‏ استفاده 
شود. مدل CA‏ مکانیزمی پردازشی است که وضعیت یک پیکسل را بر مبنای وضعیت اولیه‌اش؛ 
شرایط در پیکسل‌های پیرامونش و قوانين انتقال محاسبه می کند (داناراج و جين ۲ ) در اين بين: 
قوانین انتقال تعیین کننده‌ی رفتار سلول‌ها در زمان‌های مختلف هستند و Ol peas‏ موتور اصلی تغییرات 
در مدل CA‏ احتمال انتقال فضایی انواع کاربری زمين را تعيين می‌کنند (زنگ" و همکاران ۲۰۲۳). 
شایان ذکر است که در اين پژوهش. محدودیت‌های مکانی (مناطق ساخته شده قبلی» پهنه‌های آبی و 
پارک‌های جنگلی) و خروجی تحلیل MCDM‏ (نقشه شایستگی رشد شهری بر اساس عوامل 
محر ک) قوانین انتقال در مدل CA‏ را تشکیل می‌دهند. 

(ROC) منحنی مشخصه عملکرد نسبی‎ A Y 

ارزیابی صحت نقشه‌های شبیه‌سازی شده برای سال ۲۰۲۰ با بهره‌مندی از منحنی ROC‏ انجام 
شده است. محور عمودی اين منحنی نشانگر نرخ مثبت صحیح (TRP‏ و محور منفی اين منحنی 
نمایانگر نرخ مثبت كاذب CFRPY‏ است كه به‌ترتیب از رابطه (V)‏ و (A)‏ به دست می‌آید (فاست" 


(Ye 
_ TP 
TR د‎ A (V) رابطه‎ 
= FP 
FPR = LIN (A) رابطه‎ 


1. Isinkaralar 

2. Dhanaraj & Jain 

3. Zhang 

4. True Positive Rate 
5. False Positive Rate 
6. Fawcett 


در این پژوهش. TP‏ نشان‌دهنده تعداد سلول‌های کاربری /پوششی اراضی است که در واقعیت 
زمینی توسعه يافتهاند و مدل نيز به درستی اين گسترش را پیش‌بینی کرده است. در طرف مقابل FN‏ 
تعداد سلول‌هایی است که در حقیقت توسعه‌یافته‌انده اما مدل قادر به تشخیص آن‌ها نبوده است. FP‏ 
نيز نمایانگر تعداد سلول‌هایی است که در واقعیت توسعه نیافتند. اما مدل به اشتباه توسعه فضایی آن‌ها 
را پیش‌بینی کرده است. درنهایت TN‏ نشانگر تعداد سلول‌هایی است که در حقيقت زمینی توسعه 
نداشته‌اند و مدل نیز به‌درستی فقدان توسعه آن‌ها را پیش‌بینی کرده است. گفتنی است که مساحت 
سطح زیر اين منحنی که AUC!‏ نام دارد. بیانگر صحت مدل است و هر جه اين مقدار به عدد یک 


نزدیک‌تر باشد» صحت پیش‌بینی‌ها نیز بیشتر خواهد بود (الرفات و ليوك ۲۰۲۲). 


4 یافته‌های تحقیق 


.٤‏ ۱. اهمیت عوامل محرک و شایستگی رشد شهری 

يس از شناسایی و تعیین عوامل محرک توسعه فضایی کلان‌شهر قم ضرایب وزنی اين عوامل 
براساس نظرهای تیمی هفت‌نفره از دانشگاهیان خبره در حوزه‌های برنامه‌ریزی شهری و علوم 
جغرافیایی تعیین شده است. فرایند وزن‌دهی معیارها توسط خبرگان با روش BWM‏ انجام شده 
است که نتایج آن در جدول ۵ آمده است. گفتنی است که اين نتایج هم در فرایند اعتبارسنجی مدل و 
هم در فرایند شبیه‌سازی تغییرات LULC‏ سال ۲۰6۰ استفاده شده است. براساس نتایج وزن‌دهی؛ 
معیار فاصله از مناطق انسان‌ساخت بیشترین و معیار ارتفاع» کمترین ضریب را به خود احتصاص داده 
است. در پژوهش‌های متعددی ازجمله حسین شبهو و اسلام (۲۰۲۰) و کمراج و رنگارجان" (۲۰۲۲) 
نیز همین معیار (فاصله از مناطق انسان‌ساخت) بیشترین ضریب اهمیت را در مسئله توسعه فضايى 


شهر به خود اختصاص داده است که اين امر نشان‌دهنده اهمیت دوچندان اين شاخحص است. 


1. Area Under the Curve 
2. Al Rifat & Liu 
3. Kamaraj & Rangarajan 
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حدول 0 نتایج وزد‌دهی عوامل محر ک رشد شهری کلان‌شهر قم با روش BWM‏ 
مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


(wj) اوزان‎ 


YI dle ee‏ بط 
Bee‏ الا pe‏ 


در ادامه به کمک لایه‌های GIS‏ استانداردشده معیارها برای سال‌های ۲۰۱۰ و ۲۰۲۰ و همچنیر 


اوزان به‌دست‌آمده از نظرهای خبرگان در جدول ۵» نقشه شایستگی رشد شهری منطقه برای سال‌های 
مذکور angi‏ شده است (شکل ۵). نتایج نشان می‌دهد که سطح وسیعی از منطقه از پتانسیل زیادی 
برای توسعه فضایی مناطق شهری در آینده برخوردار است؛ بااین‌حال» مناطقی در غرب. شمال شرق و 
جنوب منطقه نیز وجود دارند که پتانسیل توسعه فضایی شهر در آنجا اندک است. توزیع فضایی 
عوامل محرک رشد شهری. گواهی بر جهت‌گیری يا فقدان جهت‌گیری احتمالی فضاى شهر در اين 
مناطق است. مقایسه اين نتایج با تحقیقی که پیش‌تر سلیمانی و همکاران (۱۳۹۹) در راستای تحلیل 
فضایی توسعه کلان‌شهر قم انجام دادند نیز نشان می‌دهد که احتمالاً جهات توسعه اين کلان‌شهر به 
همین سمت‌وسو خواهد بود. از طرفی پتانسیل توأمان رشد شهری در اين منطقه نيز نسبت بسیار 


شهرنشینی در اين کلان‌شهر )10 درصد) را توجیه می کند. 


سال بيست و يكم بهبود مدل‌سازی توسعه فضایی شهرها با تلفیق ... ۳۹ 


شایستگی رشد شهرى ۲۰۱۰ شايستكى رشد شهرى ۲۰۲۰ 


شکل ۵. نقشه‌های شایستگی رشد شهری برگرفته از تحلیل MCDM‏ 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


.٤‏ ۲. شییه‌سازی رشد شهری 
در راستای شبیه‌سازی تغییرات LULC‏ منطقه موردمطالعه برای سال ۲۰۲۰ ابتدا نقشه‌های 
پتانسیل تبدیل کاربری‌ها به‌صورت جداگانه به کمک الگوریتم‌های ANN‏ و angi SVM‏ شد. 
بدین‌منظور و برای آموزش هرکدام از این مدل‌ها از نقشه شایستگی رشد شهری حاصل از 
تحلیل MCDM‏ در سال ۲۰۱۰ استفاده شده است (شکل ۵). پارامترهای مربوط به shal‏ مدل‌های 
ANN‏ و SVM‏ در جدول + نشان داده شده است. 
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جدول I‏ پارامترهای حاصل از اجرای مدل ANN‏ و SVM‏ برای شبیه‌سازی تغییرات LULC‏ سال ۲۰۲۰ 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


ANN 


ee e 
با مس اس‎ 


در مرحله بعدی تغییرات کمی کاربری‌ها در حد فاصل زمانی بين سال‌های ۰ تا ۰ با 


استفاده از زنجیره Markov‏ محاسبه شده است (جدول AV‏ براساس اين ماتريس» 44 درصد از 


اراضی شهری (انسان‌ساخت») و ٩۰‏ درصد اراضی غیرشهری (زراعت و باغات. پارک جنگلی؛ 
پهنه‌های آبی و بایر) بدون تغییر باقی مانده‌اند. در ادامه نیز تخصیص مكانى ' کاربری‌ها با مدل CA‏ 
انجام شده و درنهایت با بهره‌گیری از هرکدام از اين الگوریتم‌های ترکیبی ANN-CA-Markov)‏ و 
۸۷ )-]۱۷ 9۱۷ دو نقشه برای LULC‏ منطقه در سال ۲۰۲۰ شبیه‌سازی شده است (شكل 


A 


جدول V‏ ماتریس تغییرات کمی کاربری‌های شهری حاصل از زنجیره Markov‏ در حد فاصل زمانی بين 
سال‌های ۲۰۰۰ تا ۲۰۱۰ 
مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


ی وهر 


1. Land Allocation 


SVM-CA-Markov ANN-CA-Markov 
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شکل 1 نقشه‌های LULC‏ شبیه‌سازی شده برای سال ۲۰۲۰ 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


: ۳. صحت‌سنحی شبیه‌سازی 

اعتبارسنجی نتایج حاصل از شبیه‌سازی تغييرات LULC‏ (شكل » با نقشه واقعیت زمینی 
LULC‏ منطقه در سال ۲۰۲۰ (شکل (E‏ با استفاده از منحنی ROC‏ انجام شده است. نتایج اين 
منحنی نشان می‌دهد که مقدار AUC‏ برای هریک از الگوریتم‌های ترکیبی ANN-CA-Markov‏ و 
SVM-CA-Markov‏ به‌ترتیب ۰/۹۵ و ۰/۹۲ است. برمبنای اين ارقام هر دو الگوریتم. توسعه 
فضایی منطقه موردمطالعه را با صحت بسیار زیادی شبیه‌سازی کرده‌اند؛ با وجود این ازآنجاکه 
الگوریتم ترکیبی SVM-CA-Markov‏ صحت بیشتری را به خود اختصاص ols‏ است. این 
الگوریتم به‌عنوان الگوریتم بهینه برای مدل‌سازی توسعه فضایی کلان‌شهر قم در چشم‌انداز ۲۰۶۰ 
برگزیده شده see‏ شايان ذكر است كه بيش تر نيز رينو و گوتزکه' (۲۰۱۵) در تحقيق خود کارایی 


در نمونه‌ای دیگر نیز فنك و همکاران »٠*(‏ کارایی این الگوریتم تركيبى را اثبات کردند؛ 


1. Rienow & 6 


alee yyy‏ جغرافیا و توسعة ناحیهای شمارة چهل و چهارم 


بااین‌حال پژوهش حاضر نشان می‌دهد که تلفیق اين مدل با زنجیره Markov‏ می‌تواند صحت انجام 


فرایند مدل‌سازی را بیش ازپیش افزایش دهد. 


ANN-CA-Markov SVM-CA-Markov 
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شکل V‏ نمودار منحنی مشخصه عملکرد نسبی (ROC)‏ 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


ء ۶. مدل‌سازی کاربری / پوشش اراضی سال ۲۰۶۰ 

پس از تعیین الگوریتم SVM-CA-Markov‏ به‌عنوان الگوریتم بهينه. از این الگوریتم برای 
شبیه‌سازی تغییرات LULC‏ در افق ۲۰۶۰ استفاده شده است. صحت و اندازه مهارت حاصل از 
اجرای الگوریتم SVM‏ به‌ترتیب مقادیر ۰/۹۱۳۹ و ۰/۸۷۵7 به دست آمد. در ادامه ماتریس تغییرات 
کمی حاصل از تشکیل زنجیره Markov‏ در جدول ۸ و نتایج تخصیص مکانی کاربری‌ها با مدل 
CA‏ و نقشه شبیه‌سازی تغییرات LULC‏ کلان‌شهر قم برای افق ۲۰۶۰ در شکل ۸ نشان داده شده 


البتت. 


جدول ۸ ماتریس تغییرات كمى کاربری‌های شهری حاصل از زنجیره Markov‏ در حد فاصل زمانی بين 
سال‌های ۲۰۱۰ تا ۲۰۲۰ 


مأخذ: یافته‌های پژوهش ۱۶۰۲ 


رت ی مرت 
مرب 


شکل A‏ نقشه LULC‏ شبیه‌سازی شده برای سال ۲۰۶۰ با الگوریتم SVM-CA-MARKOV‏ 


مأخذ: یافته‌های پژوهش, ۱۶۰۲ 


ارزیابی روند تغییرات LULC‏ سال ۰ با نتایج بيش بينى شده برای سال ۰ (جدول 4( 
نشان می‌دهد که در اثنای اين مدت توسعه فضایی کلان‌شهر قم ادامه خواهد داشت و وسعت اراضی 


شهری در اين منطقه از ۱۳۹/۸۲ کیلومترمربع در سال ۲۰۲۰ به ۱۸۳/۵۹ کیلومترمربع در سال ۲۰6۰ 


alee ۲٤‏ جغرافیا و توسعة ناحيداى شمارة چهل و چهارم 


افزايش خواهد یافت. به همین نسبت حدود ٤‏ درصد از وسعت مناطق غیرشهری منطقه نيز کاهش 
خواهد یافت. در بررسی روند توسعه فضایی کلان‌شهرهای ایران که منصوریان و همکاران (۱۶۰۰) 
انجام دادند. کلان‌شهر قم در دو دهه اخیر با ۳/۲ درصد بعد از کلان‌شهرهای شیراز و تبریز» بیشترین 
نرخ رشد فیزیکی را به خود اختصاص داده است. همچنین تجزیهوتحلیل فضایی روند توسعه 
کلان‌شهر قم نشان داد که توسعه فضایی اين شهر غالباً به سمت جنوب و جنوب غرب منطقه بوده كه 
اين بخش base‏ شامل اراضی زراعی-باغی و نواحی روستایی است. روستاهای براسون قلعه چم و 
زرقان ازجمله مناطق روستایی هستند که در اين بخش واقع شده‌اند. توسعه فضایی توأمان منطقه به 
سمت اراضی زراعی و نواحی روستایی می‌تواند آسیب‌های زیست‌محیطی و اقتصادی زیادی را به اين 
منطقه وارد کند؛ چراکه اين امر از یک سو باعث کاهش سطح زیرکشت محصولات زراعی شده و از 
سوی دیگر پیوستن مناطق روستایی به پیکره شهر, می‌تواند منجر به کاهش سرانه فضای سبز و 
امکانات رفاهی در اين کلان‌شهر شود. در اين راستاء اتخاذ استراتژی‌های هدفمند و يايش جهات 


توسعه شهر از سوی سازمان‌ها و متولیان مربوط بسیار ضروری به نظر می‌رسد. 


جدول 4. مقایسه تغييرات LULC‏ نقشه واقعیت زمین ۲۰۲۰ با نقشه شبیه‌سازی شده سال ۲۰۶۰ 


مأخذ: یافته‌های پژوهش, ۱۶۰۲ 
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(KM2) 14 ۱۳۹ 
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سال بيست و يكم بهبود مدل‌سازی توسعه فضایی شهرها با تلفیق ... ۳۳۹ 
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کارت اد دوس 


شکل 4 رشد شهری کلان‌شهر قم از سال ۰ تا سال ۲۰۶۰ 


مأخذ: یافته‌های پژوهش. ۱۶۰۲ 


0. نتيجه گیری و پيشنهادها 


یکی از اصلی‌ترین ضروریات در تصميم كيرىهاى مربوط به امور شهری. داشتن اطلاعات به 
هنگام از تغییرات کاربری/پوشش اراضی است. مسلم است. در اختیار داشتن چشم‌انداز توسعه فضایی 
شهرها متولیان مربوط را قادر می‌سازد که زیرساخت‌های لازم را فراهم آورند. در اين راستا و با هدف 
ارائه چشم‌اندازی جامع از توسعه فضایی کلان‌شهر قم در افق ۲۰۶۰ پژوهش حاضر به مدل‌سازی 
توسعه فضایی اين کلان‌شهر و حریم پیرامون OF‏ پرداخته است. در اين راستا و به‌منظور دستیابی به 
الگوریتم بهینه برای بهبود عملکرد مدل‌سازی, دو الگوریتم یادگیری ماشین پرکاربرد در زمینه 
موردبحث یعنی ANN‏ و SVM‏ با مدل ترکیبی Gali CA-Markov‏ شده و تغییرات LULC‏ 
کلان‌شهر قم برای سال ۲۰۲۰ شبیه‌سازی شده است. نتایج نشان داد که الگوریتم ترکییی SVM-‏ 
CA-Markov‏ با مقدار AUC‏ معادل ۰/۹٦‏ از صحت و تعمیم‌پذیری بیشتری در مقایسه با الگوریتم 


alee ۳۳۱‏ جغرافیا و توسعة ناحیه‌ای شمارة جهل و چهارم 


ANN-CA-Markov‏ برخوردار است. همین امر سبب شده است که مدل‌سازی توسعه فضایی 
منطقه در افق ۲۰۶۰ با استفاده از اين الگوریتم انجام شود. گفتنی است که نقشه‌های چندزمانه 
LULC‏ و نقشه شايستكى رشد شهری برگرفته از تحلیل MCDM‏ ورودی‌های مدل را تشکیل 
می‌دهند. نتایج حاصل از اجراى مدل بركزيده نشان داد كه وسعت اراضى شهرى کلان‌شهر قم تا سال 
۰ بیش از ۳۱ درصد افزايش خواهد داشت و به بيش از ۱۳۸ كيلومتر مربع خواهد رسيد. 

نتایج اين پژوهش می‌تواند پنجره‌ای رو به آينده برای مدیران و متولیان امور شهری کلان‌شهر قم 
بگشاید و ایشان را در راستای فراهم آوردن زیرساخت‌های لازم. تدوین سیاست‌های منطقه‌ای در 
زمینه توزیع متوازن جمعیت و خدمات. برنامه‌ریزی و هدایت سرمایه گذاری و نيل به سمت تحقق 
اهداف توسعه پایدار يارى Whey‏ در انتها پيشنهاد می‌شود که در مطالعات آتی» ضمن بهره‌مندی از 
معیارهای تصمیم‌گیری جامع‌تر نتایج الگوریتم مورداستفاده در اين پژوهش نيز با سایر الگوریتم‌های 


مدل‌سازی رشد شهری مقایسه شود تا دقت نتایج بیش ازپیش افزایش یابد. 
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